Generation of broad XUV continuous high harmonic spectra and isolated attosecond 

pulses with intense mid-infrared lasers 
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We present experimental results showing the appearance of a near-continuum in the high-order 
harmonic generation (HHG) spectra of atomic and molecular species as the driving laser intensity 
of an infrared pulse increases. Detailed macroscopic simulations reveal that these near-continuum 
spectra are capable of producing lAPs in the far field if a proper spatial filter is applied. Further, our 
simulations show that the near-continuum spectra and the lAPs are a product of strong temporal 
and spatial reshaping (blue shift and defocusing) of the driving field. This offers a possibility of 
producing lAPs with a broad range of photon energy, including plateau harmonics, by mid-lR laser 
pulses even without carrier-envelope phase stabilization. 



It has been a decade since the first reported genera- 
tion of isolated attosecond pulses (lAPs) l] and the first 
full characterization of an attosecond pulse train (APT) 
0, both seminal results that gave birth to the field of 
attosecond science 0-13. While APTs can be now rou- 
tinely generated based on high-order harmonic genera- 
tion (HHG) in gases, production of lAPs is still a very 
active field of research. Common to all methods of gen- 
erating lAPs is the concept of temporal gating, which is 
used to select a XUV burst from an APT. Current meth- 
ods can be divided into three different categories. The 
most intuitive, but technically demanding, approach is 
based on XUV spectral filtering near the cutoff of the 
HHG spectra. This approach requires a short driving 
pulse (about two optical cycle) with a stabilized carrier- 
envelope phase, d, [I]. This approach could be called 
intensity gating, as it selects the harmonics emitted near 
the peak of the most intense half-cycle of the (short) 
laser pulse. The second approach is based on polariza- 
tion gating [gj-fllj , and the generalized double optical gat- 
ing |12h14| , in which the short driving pulse requirement 
is somewhat relaxed. While these two approaches have 
their origin already in the single-atom response, the third 
approach is based on the macroscopic propagation ef- 
fects |25i], which includes, ionization gating [16- 18], phase 
matching [Toj , and ionization-driven reshaping Gaarde et 
al [H, [131 ■ -"-t ti^s ^1^0 been demonstrated recently that 
near-continuum spectra can be generated with two-color 
multicycle driving laser pulses by carefully adjusting the 
wavelength of the supplementary pulse ^] . 

In order to produce a short lAP, XUV radiation with a 
broad frequency bandwidth needs to be generated. Since 
the HHG cutoff law is ilcutof / ^ there are only two 

ways to extend Qcutof f '■ increasing the laser intensity II , 
or increasing the laser wavelength A l ■ The first option is 
limited by ionization, which dramatically reduces HHG 
yield due to the depletion of neutral atoms and, more 



importantly, the phase mismatching due to free electrons 
in the medium. Using a short laser pulse does not help 
much, since the current technology nearly reaches the 
one-cycle limit already. The second option is to use lasers 
with a long wavelength. This comes with an additional 
advantage of generating a shorter lAP since the atto- 
chirp is inversely proportional to the wavelength of the 
driving laser The price to pay is that HHG yields 

have been shown to drastically decrease with wavelength 
[23 - |25| . However, it has been demonstrated recently that 
this decreasing yield with laser wavelength at the single 
atom level can be compensated for by increasing the den- 
sity of atoms resulting in increased XUV flux for some 
phase matching conditions [26{ . 

In this paper we present experimental evidence of gen- 
eration of a broad near-continuum in HHG spectra from 
atomic and molecular gases by using two-cycle mid-IR 
(1825 nm) laser pulses from atomic and molecular tar- 
gets. As supercontinuum emission is consistent with 
an lAP, we carry out detailed theoretical simulations to 
show that this scheme indeed is capable of producing 
lAPs if a proper spatial filter is applied. Our analysis 
shows that the generation of a near-continuum and an 
lAP is the result of combined effect of reshaping of the 
driving laser and non-adiabatic phase mismatch, caused 
by the free electrons due to sub-cycle strong field ioniza- 
tion. With this new method and the light source used, it 
is possible to use a broad range of harmonics, including 
plateau harmonics, to generate an lAP with a tunable 
XUV photon energy. 

The experimental setup as shown in Fig. [1] starts with 
a Ti:Sa 800nm, 35 femtosecond laser delivering 6 mJ of 
energy per pulse, at a repetition rate of lOOHz. Using 
the idler beam of a fiuorescence-secdcd commercial opti- 
cal parametric amplifier (OPA, HE TOPAS from Light 
Conversion) the 800nm pulses are converted to 1825 nm. 
The idler signal from the OPA produces pulses with 52 



100 HzTi:Sa 



OPA 



950 laJ 



1.8 |j,m 



R=2m 



hollow-core fiber f= Im 



FS P 



gas jet MC 
\J ■ 




R= 0.5 m slit grating 



FIG. 1: Experimental setup. 



fs duration and 900 /i J which are then focused into a 400 
/im hollow core fiber. Spectral broadening occurs in the 
capillary due to self-phase modulation and the pulses are 
compressed by linear propagation through bulk material 
(27| . The pulse duration, measured with a second har- 
monic generation frequency resolved optical gating device 
yields 14 fs, corresponding roughly to a 2 cycle pulse. 
Following the fiber, the beam is collimated with a Im 
lens and then focused in a thin (500 [im) pulsed jet of 
atoms or molecules. Before entering the chamber the 
beam passes through a half wave plate and a nanoparti- 
cle thin film polarizer (Thorlabs LPMIR) that allows for 
intensity control while keeping all other optical param- 
eters constant. Intensity calibration is done using the 
method outlined in p3| . In short, we measure the ions 
in parallel with the harmonic emission as a function of 
laser energy and then fitted to a space and time aver- 
aged Yudin-Ivanov ionization model 28] in a cylindrical 
geometry [2^ . The resulting equations are only intensity 
dependent once the pulse duration and wavelength are 
known. All other pulse and detector parameters enter as 
a single scaling coefficient. Our XUV spectrometer uses 
a 1200 lines/mm ruled grating (Hitachi), which disperses 
the light into a microchannel plate/phosphor screen de- 
tector (Burle APD 3115 32/25/8 I EDR MgF2 P20). An 
estimate of the spectrometer resolution can be made by 
looking at the recombination lines of ionized O2 using 
circularly polarized light [sO]. Through this method we 
estimate an spectral resolution of < \eV in the range of 
30eV to 70eV. 

In Fig. [5] we show the HHG spectra obtained in atomic 
Xe and molecular NO at different laser intensities. The 
most striking feature is that the harmonic spectra be- 
comes a near-continuum as laser intensity increases. This 
effect is particularly strong for low photon energies (51 < 
40 eV). These continuous spectra extend over a broad 
range of photon energy from the cutoff at about 100 eV 
down to very low energy of 20 - 30 eV. We also observe 
this feature in other atomic targets as well as in other 
molecules such as N2, O2, CO, CO2, and other organic 
molecules with relatively low ionization potentials. We 




40 60 80 
Photon energy (eV) 



100 



FIG. 2: Experimentally measured HHG spectra in Xe (a) 
and NO (b) for dilferent laser intensities. Driving field has 
\l = 1825nm and duration of 14 fs. Intensity values are in 
units of Jo = lO^T^/cm^ 



therefore believe that this behavior is quite universal and 
is independent of atomic and/or molecule species. We 
also note a saturation effect in observed HHG spectra, 
when further increase of laser intensity does not change 
HHG spectrum significantly. This can be clearly seen 
in case of Xe for laser intensities above about 2 x 10^'* 
W/cm^ in Fig. 1(a). 

As the appearance of a broad continuum is an indica- 
tion of a possible lAP, it is very tempting to speculate 
this transition as a transition from an APT to lAP by 
controlling the intensity of few-cycle mid-IR pulses. The 
transition from APT to lAP at very high laser intensities 
has has been demonstrated previously where the mecha- 
nism was ionization gating [13, Hlj . In ionization gating, 
the field is so intense that the electron is momentarily 
stripped from the atomic core, thus making the electric 
field effectively one-cycle. To confirm our hypothesis of 
the transition from APT to lAP and to further under- 
stand the underlying physics we carried out theoretical 
simulations, which include the macroscopic propagation 
of both driving laser and HHG fields, see below. Our 
simulations reveal that the mechanism behind this tran- 
sition is quite different from the conventional ionization 
gating. In fact, the mechanism is closer to the ionization- 
driven reshaping of the driving laser pulse, as discussed 
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by Gaarde et al [l5|,l20|. However, the use of mid-IR laser 
pulse leads to a different regime, where the gas pressure is 
much lower, and the gas Jet is much thinner, as compared 
to that of Gaarde et al [3) HOl ■ 

Before presenting results of our simulation we first 
comment on the differences between the experimental 
conditions for the present work and for our earlier re- 
sults [S^l- The conditions for 32 1 were chosen to avoid 
phase mismatch so as to measure harmonic spectra that 
come as close as possible to the single atom response. 
In doing so a broad enhancement in HHG yield was ob- 
served near 100 eV, which is attributed to the effect of 
the interchannel coupling in the photoioinzatin / photo- 
recombination cross section of Xe ^32]. To access the 
single atom response the laser was focused l-2mm down- 
stream from the jet where the gas density was relatively 
low supporting harmonics which were well phase matched 
across the entire xuv spectrum. In the present work, the 
focus was brought closer to the jet into a region where 
the gas density was higher. This lead to the sub-cycle 
modification of the laser field that we discuss in the fol- 
lowing section and also explains why the cutoff observed 
after macroscopic propagation falls short of the cutoff ex- 
pected based on the laser intensity and the single atom 
response. The laser setup used for the present work is 
identical to that in [s^ although the pulse duration is 
somewhat longer (14 fs vs llfs). The small bump near 
40eV for the highest intensity of 2.0 x 10^* W/cm^ in 
FigHh-) is an artifact of the grating used for the present 
study. To explain the spectrum that we observe in the 
present work we performed simulations in atomic Xe at 
different focusing positions, gas pressure, and pulse du- 
rations. We found that the most influential parameters 
are the pressure and the size of the gas jet 

For an in-depth detail of the theoretical method used 
in this paper, we point the reader towards Ref. [33| . 
We only briefly describe the method here. The theoret- 
ical study of high harmonic generation (HHG) consists 
of two parts _34-36]. The first one is to calculate the in- 
duced dipole of each atom or molecule in the laser field. 
The second part is to consider the nonlinear propagation 
of the fundamental laser pulse and the harmonic fields 
in the medium by solving the Maxwell's equations. For 
the propagation of the fundamental laser field, we in- 
clude the effects of refraction, nonlinear Kerr effect, ion- 
ization, and plasma defocusing [3^. Due to the absence 
of such parameters for most molecules, our simulations 
focus mainly on atomic Xe. We use the recently devel- 
oped quantitative rescattering (QRS) theory |37H39j to 
calculate the induced dipole of single atom or molecule. 
Once the induced dipoles are calculated for different laser 
intensities in the medium, they are then fed into the prop- 
agation equations of harmonic field. At the exit face of 
the medium, we obtain the near-field harmonic emission. 
Harmonics then propagate further in the free space, and 
they may go through a slit before they are collected by 
the detector. The far-field harmonic emission can be ob- 
tained from the near-field harmonic emissions through 
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FIG. 3: (color online) HHG spectra calculated using QRS the- 
ory for different peak intensities as indicated where /o^lO^* 
W/cm^. Panel (a) shows spectra for a single atom response. 
Panel (b) shows CEP-averaged spectra from macroscopic sim- 
ulations for different intensities. Panel (c) shows spectra from 
macroscopic simulations for gas pressures of 10 and 25 Torr, 
at a CEP of n/2. Other laser parameters are the same as in 
the experiment for all panels. 



a Hankel transformation. We pay particular attention 
to the multi-electron correlated process that occurs in 
atomic Xe for returning electrons with energies beyond 
90 eV. The QRS states that the single-atom response 
induced dipole D{uj) can be written as 



D{uj) = W{cu)diLo), 



(1) 



with W{u}) the laser-driven macroscopic wave packet and 
d{Ld) the photo-recombination (PR) dipole moment. To 
describe the PR process of Xe, we need to use the PR 
transition dipole moment including multielectron effects 
in QRS theory. Such PR moment can be extracted from 
the experimental photoionization cross section (PICS) 
[40I . I4H . In the calculation, we incorporate the PICS 
of Xe, which is extended into high-photon energy region 
(a few hundred eV), using the relativistic random-phase 
approximation (RRPA) by Kutzner et al. . Note that 
the phase of PR transition dipole moment is assumed the 
same as that solely caused by 5p. 
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The results of our simulations in atomic Xe are shown 
in Fig. [31 First, the CEP-averaged HHG spectra of the 
single atom response obtained with the QRS method are 
shown in Fig. Ela) for two laser intensities of 0.5 x 10^'' 
W/cm^ and 1.0 x lO^'' W/cm^. Clearly, the simulation 
at the single atom level does not reproduce the the tran- 
sition to the continuum spectrum, observed in the exper- 
iment, as well as the HHG cutoff position for the high 
laser intensity. The ionization probability is about 10- 
30% at the end of the pulse in this range of laser in- 
tensity (Jl « 1.0 X 10" W/cm^). This result basically 
excludes the possibility of ionization gating at the sin- 
gle atom level, as a mechanism for the transition to the 
near-continuum spectra. We also found that further in- 
crease in laser intensity up to 2.0 x lO^'' W/cm^ does 
not change this behaviour. We comment that the genera- 
tion of an lAP with a high pulse energy (few nanojoules) 
has been reported recently by Ferrari et al fsTj, which is 
based on complete depletion of the neutral Xe atoms. 

Next, we show in Fig. [2Ib) the harmonic spectra for 
four peak laser intensities after propagation through the 
gas jet with parameters that closely resemble the exper- 
imental conditions. Such condition include not only the 
laser parameters but also the jet size (1 mm), the spec- 
trometer slit opening (190 /im) and distance of the slit to 
the interacting region (455 mm). In the theoretical calcu- 
lations the beam is assumed to be Gaussian. All spectra 
are averaged over random values of the CEP as in the 
experiment. The macroscopic simulations indeed show 
a transition to a near-continuum HHG spectrum for a 
broad range of photon energy from the "apparent" cutoff 
down to about 30 eV, as the laser intensity increases, thus 
reproducing nicely the experimental finding. The simu- 
lations also reproduce the saturation effect observed in 
the experiments. This is not the case for the single atom 
simulations shown in Fig. |31[a) . This lead us to conclude 
that the near-continuum is a consequence of the propa- 
gation in the medium of both the fundamental and the 
harmonics fields. 

By comparing HHG spectra obtained for the single 
atom with the ones from a macroscopic simulation for the 
same peak intensities, we observe large discrepancies in 
the cutoff and the overall shape. Before explaining these 
discrepancies as well as the physical mechanism behind 
the experimentally observed spectra we address here the 
differences between the current experimental data and 
that of Shiner et al [32|. Fig. 3(c) shows a comparison 
between the HHG spectra from macroscopic simulations 
for two different gas pressures of 10 and 20 Torr. The 
laser intensity is 2.0 x 10^"* W/cm^ and the CEP is cho- 
sen to be 7r/2. The "apparent" cut-offs are different by 
about 15 eV. Thus the discrepancy between the two ex- 
periments can be tentatively attributed as due to the 
different gas pressures. Moreover, we note that the size 
of the gas jet could also affect the cut-off position and the 
overall shape of HHG spectra. A thinner gas jet would 
give a spectrum closer to a single-atom response, i.e., 
with a higher cut-off. 
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FIG. 4: Driving electric field (Al = 1825nm) as a function of 
time (in units of the optical cycle) at three different positions 
relative to the atomic jet, for a gas pressure of 25 Torr (a) 
and 10 Torr (b). All values are "measured" on axis. Laser 
parameters are II = 2.0 x 10^'' W/cm^ and CEP=7r/2. 



In order to understand the mechanism behind the tran- 
sition to near-continuum HHG spectrum, we followed the 
fundamental field in space and time through the simula- 
tion. In Fig. |4] we show the on-axis electric field ampli- 
tude as a function of time at the entrance, middle, and 
just before the exit of the Xe gas jet, for a gas pressure of 
25 Torr and 10 Torr. The CEP is chosen to be 7r/2, and 
the laser intensity is 2.0 x 10^^ W/cm^. As can be seen 
from the figure, the temporal shape of the laser pulse is 
strongly modified as it propagates though the medium. 
In particular, a dynamic blueshift (or chirp) can also be 
seen clearly. In fact, this shape is quite typical for the 
self-phase modulation found in a rapidly ionizing medium 
[isl [20I [33I I . A chirped pulse means that harmonics will 
be produced at different fundamental frequencies depend- 
ing on where and when they are generated. The effect 
of chirp in HHG has already been established [i^ HJ] 
and it induces a frequency shift in the observed harmonic 
spectrum. In our case, the harmonics are produced with 
different chirped pulses all at once. Furthermore, the 
electron density strongly varies as a function of radial 
distance due to the radial variation of the laser intensity. 
This radial variation electron density acts like a negative 
lens to defocus the laser beam. This has also been found 
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in the simulation (see, [s^). We found that both tem- 
poral and spatial reshapings are quite negligible at the 
laser intensity of 0.5 x 10^"* W/cm^, since the ionization 
is insignificant. This explains the observed transition to 
the near-continuum in HHG spectra as laser intensity in- 
creases. 

The temporal and spatial reshaping of the laser beam 
significantly reduces the laser intensity in the medium, 
which accounts for the "apparent" lower energy cut-off 
in the HHG spectra found in the macroscopic simula- 
tions, as compared to the single atom simulations. Fur- 
ther increase in laser intensity above 2.0 x 10^* W/cm^ 
will leads to even stronger reshaping and redistribution 
of laser intensity, which explains the saturation effect dis- 
cussed above. 

The reshaping of the laser pulse has been discussed 
before in the context of the 800 nm driving pulse f5, 20]. 
The effect is even more significant for the mid-IR laser, 
since the effective index of refraction due to free electron 
is ~ —n^Xj^, where is the density of free electron 
and Al is the wavelength of the driving laser. Therefore 
the reshaping of the laser pulse occurs at much lower gas 
pressure and size (25 Torr, 1 mm long) as compared to 
that of Gaarde et al (135 Torr, 3 mm long). We further 
comment that at the lower pressure of 10 Torr, shown 
in Fig. 4(b), the reshaping is less severe as compared to 
25 Torr, thus leads to the higher HHG cutoff seen in 
Fig. 3(c). 

To further understand the nature of the near- 
continuum and to address the usefulness of our scheme for 
the production of lAPs we follow the attosecond pulses in 
time as they emerge from the interacting region. We use 
a spectral filter to select harmonics from H40 (27.2eV) to 
H80 (54.4 eV). Fig.[5l[a) shows the attosecond pulses us- 
ing these harmonics for the same driving field [Xl = 1825 
nm and duration (FWHM) of 14 fs] with a peak inten- 
sity of 1.5 X lO^"* W/cm2, and a CEP of 7r/2 at the end 
of the interacting region (near field). The presence of 
an APT is evident. As proposed in [15|, an lAP can be 
obtained in the far field by spatially filtering the XUV 
pulses. In fact, due to the temporal and spatial reshap- 
ing, the HHG emission at different times will have a dif- 
ferent divergence. Fig. [Sjb) shows harmonics H40 to H80 
in time domain, under the same conditions as panel (a), 
but the pulse is "measured" in the far field after a circular 
filter with a radius of 300 /im, located 455mm after the 
focus. In this case, a single 390 as pulse is indeed mea- 
sured, thus proving that the observed XUV supercontin- 
uum can indeed produce lAPs. This attosecond burst 
is generated in the leading edge of the pulse, less than 
one optical cycle before the peak intensity is reached. 
Clearly, the more divergent attosecond bursts, generated 
at later times have been blocked by the filter. This result 
also indicates that the phase mismatching due to the free 
electron might also play an important role, as the near- 
axis contribution from the later sub-cycles is suppressed 
due to the phase mismatch. We mention that the phase 
mismatch due to the free electron is more severe for the 
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FIG. 5: (Color online) (a) Intensity of attosecond pulses in the 
near field, (b) Intensity of isolated attosecond pulse (lAP) in 
the far field (measured 455 mm after the laser focus). Black 
is using a spatial filter with a radius ro = 300/im, located 
455mm after the focus. In green, a spatial filter of To = 
3mm, at the same position is used instead. Harmonics used 
to synthesize attosecond pulses are II40 (27.2 eV) to H80 (54.4 
eV). Laser intensity at the focus is 1.5 x 10^* W/cm^ and the 
CEP is 7r/2. 



long wavelength. Figure Eljb) also shows in green (color 
online) the XUV pulse under the same conditions but 
using a 3mm spatial filter. In this case, an lAP is not 
obtained, thus showing the importance of the filtering. 
We also carried out simulations with a lower laser inten- 
sity of 0.5 X 10^^ W/cm^. The result shows an APT even 
after the spatial filter is applied. At this low intensity, 
the ionization and the pulse reshaping are insignificant. 
Therefore, in this case different XUV bursts at different 
half-cycles will have similar divergence. 

We emphasize that relatively low harmonics is chosen 
here to illustrate the usefulness of our scheme for pro- 
ducing an lAP with plateau harmonics. In contrast, in 
a typical lAP generation with the reshaping mechanism 
with 800 nm laser pulses, only the harmonics near the 
cutoff are used for producing an lAP [T^. In Ref. [3^ 
we provide further details on how this lAP can be ma- 
nipulated by changing the CEP of the driving field, as 
well as the spectral and spatial filters. In particular, it 
was shown that the lAPs can be produced even with- 
out CEP stabilized laser pulses for most of the CEPs. 
Although the experimental results and theoretical analy- 
sis presented here are for two-cycle pulses, we anticipate 
that production of lAPs within our scheme is also feasible 
with longer pulses. 

In conclusion, we have observed a near continuum in 
the HHG spectrum, with a 1.8 /im laser driver as the 
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laser intensity increases. By performing theoretical sim- 
ulations of the full macroscopic HHG process combined 
with a realistic single atom response calculation we are 
able to reproduce the same behavior as we observe in the 
experiment. The simulations show that the transition 
from a train of attosecond pulses to an isolated attosec- 
ond pulse is due to strong reshaping of the driving laser as 
it propagates through the medium. The main advantage 
of our approach for producing lAPs is that the central 
frequency of the XUV pulse can be tuned over a broader 
range of photon energy as compared to the standard 800 
nm. This is due to the capability of using a large range of 
harmonics, including plateau harmonics, of a broad near- 
continuum spectrum. We comment that our method is 
based on a different mechanism from that of Ferrari et al 
[Slj , which is based on complete depletion of the neutral 



atoms. However, in both methods, both plateau and cut- 
off harmonics can be used t o sy nthesize lAPs. Finally, it 
has been reported recently [2a| that the XUV yield can 
increase with wavelength when a proper phase matching 
conditions are exploited, thus creating a great venue for 
our approach. 
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